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In  this  paper  we  discuss  the  linear  regulator  problem 
associated  with  delay-differential  equations.  Formulas  for 
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in  the  state  space  pf  the  delay  equation.  Finally,  the 
linear  regulator  problem  for  the  delay  equation  is  approximated 
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In  this  paper  we  discuss  the  linear  regulator  problem  associated  with 
delay-differential  equations.  Formulas  for  the  optimal  feedback  control  and 
optimal  trajectory  are  derived;  this  naturally  leads  to  a  Riccati  integral 
equation  in  the  state  space  of  the  delay  equation.  Finally,  the  linear  regulator 
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Introduction  and  Notation. 

We  study  the  linear-quadratic  optimal  control  problem  for  a  certain  class 
of  delay-differential  equations.  In  Section  2  an  analytical  solution  in  terms 
of  a  Riccati  integral  equation  in  the  state  space  of  the  delay  equation  is 
presented.  Subsequently,  in  the  third  section,  we  develop  a  general  scheme  for 
the  approximation  of  this  optimal  control  problem.  The  averaging  approximation 
scheme  and  the  spline  schemes  are  seen  to  be  special  cases. 

As  usual  IRn  denotes  the  n-dimensional  Euclidean  space  with  norm  | * |  and 


inner  product  ( • , * ) 


IR 


n 


For  the  interval  (a,b)  c  (-“,<»)  the  Hilbert  space  of 

rb 


n  f°  2 

equivalence  classes  of  measurable  functions  x:(a,b)  -*]R  with  |x(s)|  ds  <  00 

2  n  3 

is  denoted  by  l/(a,b;JRn)  and  is  endowed  with  the  usual  inner  product  and  norm. 

2  n  2 

For  L  (-r,0-IR  ),  r  >  0,  we  simply  write  L  .  The  state  space  for  the  presenta¬ 
tion  will  be  Z  ■  lRn  *  L2  which  in  an  obvious  way  becomes  a  Hilbert  space  with 
inner  product  (•,•)  .  Generally,  norms  of  elements  will  be  denoted  by  |*| 

z 

and  operator  norms  by  ||  •  ||  .  The  set  of  all  bounded  linear  operators  between 
Hilbert  spaces  X  and  Y  is  j^(X,Y)  and  A*  denotes  the  Hilbert  3pace  adjoint 
of  the  operator  A.  Finally,  for  x:[-r,a)  -►Jl”,  a  >  0,  the  symbol  xt,  0  <_  t  <  a, 
stands  for  *t(s)  *  x(t+s)  for  sG[-r,0]. 


2.  The  Optimal  Control  Problem. 

We  are  concerned  with  functional  differential  equations  of  the  type 


ft  r° 

x(t)  -  1  A.x(t-r.)  +  A  .  (s)x(t+s)ds,  for  t  0 

i-0  1  1  J-r 

(x(0),x0)  -  (n,<p)  €  z 


(2.1) 


» 


2. 


where  A^,  i=0 . £  are  real  n  xn  matrices,  A_^  is  an  n  *  n-matrix  of 

2 

L  -valued  functions  and  0  «  r^  <  ...  <  r£  *  r.  Although  we  restrict  ourselves 

to  the  autonomous  case  here,  the  results  are  remain  true  In  the  case  where 

A  depends  on  t.  For  convenience  we  introduce  the  operator  L  given  by 
1  %  jO 

Ltp=  £  A.<P(-r.)  +  A  1(s)U(s)ds.  It  is  quite  well  known  that  (global) 
i=0  J-r 

solutions  x  =  x(*;n,4$  of  (2.1)  exist  and  that  they  do  not  depend  on  the 

2 

representation  of  an  equivalence  class  £  L  [7];  here  x(*;n,<P)  is  called 
the  solution  of  (2.1)  if  x(s;n,49  ■  <P(s)  almost  everywhere  on  [~r,0]  and 
x(t;n,*P)  =  n  +  f  Lx  (sn,<P)ds  for  t  £  0.  By  T(t):Z  -*•  Z  we  denote  the 

Jo  8 

family  of  solution  operators  associated  with  (2.1)  via  T(t)(n,tp)  “ 

(x(t;n,49 ,xt(- ;n,4)) .  It  is  well  known  that  T(t)  is  a  linear  (^-semigroup  of 
bounded  linear  operators  for  which  ||T(t)  ||  £  M  exp(tot)  holds  for  some  M 
and  0)  in  IR;  see,  e.g.  ,[1] . 

We  can  now  state  the  following  optimal  control  problem  associated  with 

(2.1): 


|Find  u  £  L^(0,t*jRm)  which  minimizes 


ft*  ft* 

J_(n,<P,u)  -  (Fx(t*),x(t*))  +  (Dx(t) ,x(t))  dt  +  (Cu(t) ,u(t) )  dt 

0  ]Rn  Jo  Rn  Jo 

(P)  < subject  to 

x(t)  *  Lxt  +  Bu(t)  for  t  £  [0,t*],  and 

(x(0)  ,Xq)  -  (n,<P),  where  t*  >  0  and  (n,<P)  €  Z  are  given. 


The  n  x  n-matrices  F  and  D,  the  m  x  m-matrlx  C  and  the  n  x  m-matrix  B, 
which  are  assumed  to  be  independent  of  time  for  the  sake  of  a  simpler  presentation, 
satisfy:  F,D  and  C  and  are  symmetric  and  F  >  0,  D  >  0,  C  >  0.  Of  course,  (P) 


» 


3. 


is  an  extension  to  delay-differential  equations  of  the  unrestricted  finite 
terminal  time  linear  regulator  problem,  which  is  well  known  for  ordinary 
differential  equations.  But,  as  a  consequence  of  the  fact  that  the  statespace 
of  delay-differential  equations  is  a  function  space  over  the  delay  interval 
[-r,0]  (in  our  case  Z) ,  the  analysis  of  the  problem  becomes  much  more 
difficult  than  for  ordinary  differential  equations.  Many  an  earlier  inves¬ 
tigation  was  directed  towards  studying  (P)  and  we  shall  mention  some  of  them 
further  below.  Here  we  first  reformulate  (P)  in  the  state  space  Z. 


rFind  u  E  L^(0,t*;lRm)  which  minimizes 

ft* 

J0(n»<P»u)  =  <^S(t*),z(t*))z  +  (^z(t)  ,z(t))zdt  + 


rt* 

(Cu(t) ,u(t))  dt 
JO  JRm 


(^K 


subject  to 

z(t)  -  T(t)(n,<p)  + 

for  given  t*  >  0 


[  T(t-s)i?u(s)ds 
J0 

and  (n,<p)  €  Z. 


for 


t  e  [0,t*] 


(2.2) 


We  used  the  notation  5*(0,4O  “  (F0,O),  ■  (D0,O)  and  _#(6,'l0  * 

(B0,0).  In  [l]  it  is  proved  that  the  solution  z(t)  ■  z(t;u)  is  related  to 
the  solutions  x(t;n,<P,u)  of 

x(t)  -  lxt  +  Bu(t) 

(x(0)  ,xQ)  »  (T>,q>) 

via  z(t;u)  -  (x(t;n,<p,u) ,  xt(,;0»«P>u))  for  all  (n,<P)  €  Z.  Therefore,  (J?) 
and  (P)  are  equivalent  in  the  sense  that  u*  with  associated  trajectory 
x*(*;n»<Au*)  a  solution  of  (P)  if  and  only  if  u*  with  associated 
trajectory  t  (x(t;n.  <ftu*),  xfc(' ;n,  is  a  solution  of  (£*) .  The  usual 


solution  of  (P)  for  the  ordinary  differential  equation  case  (or  even  in  the 

case  of  functional  differential  equations)  involves  a  Riccati  differential 

equation  of  one  appearance  or  another  [3,  pp.130;  4,  13].  In  our  approach 

the  place  of  the  Riccati  differential  equation  will  be  taken  by  an  appropriate 

Riccati  integral  equation,  which  we  shall  explain  next.  First,  note  that  (<?) 

obviously  has  a  unique  solution  u*;  this  is  a  consequence  of  the  fact  that 
2  + 

J0(n,<P,  ’)  :L  -*-lR  is  a  strictly  convex,  continuous  and  radially  unbounded 
functional  [9] .  The  optimal  control  u*  can  be  found  from  the  equation 

Jg(n,<P,u)v  -  0  for  all  v  €  L2(0,t*-lRm) ,  (2.3) 

where  J^n.tp.u  )v  denotes  the  Frechet  derivative  at  u  in  direction  of  v. 

After  some  calculation  one  finds  that  u*  is  given  by 

u*(t)  -  -((VQ)-1W0(ri,<p))(t)  for  almost  every  t  €  [0,t*], 

where  VQ  6  ^(L^O.t*-®®)  ,L2(0,t*-Rm))  and  WQ  €  ^’(Z,L2(0,t*jmm))  are 
given  by 

VQ  -  C 

and 

wQ  -  a*f0*@T0  + 


1 


9 


5. 


Here  ^  £  ML2(0,t*;Z),L2<0,t*;Z)),  S*j  £  ^(L2(0,t*;Z)  ,Z)  and 
Tq  £  isf(Z,L2(0,t*;Z))  are  defined  by 


(J^ip)(s)  =  |  T(s-n)«p(n)dn  for  «p €  L2(0,t*;Z) 


(^0<P)  =  (^Q<P)(t*),  (TQz)(t)  =  T(t)z  for  z  €  Z, 


ft* 

C^*<P)(t)  -  j  T*(n-t)«p(n)dn,  (^0*z)(t)  =  T*(t*-t)z. 


Now  for  s  €  [0,t*]  let  us  consider  the  optimal  control  problem  (P)  ,  which 
is  defined  by  replacing  'O'  by  's'  in  (P).  Then,  similar  to  the  way  we 
found  u*,  one  can  derive  the  optimal  control  (u*)  of  (P)  which  is 

given  by  (u*)g(t)  =  ((Vg)  H?g  (n/P))(t)  almost  everywhere  in  [s,t*]  with 
Vg  €  i^L2(s,t*jlRm),  L2 ( s , t jlRm) )  and  Wg  €  .^(Z,L2(s,t*;lR,n))  defined 
analogously  to  and  W^.  For  (n,<p)  €  Z  the  optimal  trajectory  corres¬ 

ponding  to  (P)g  is  therefore  given  by 


S(t,s)(n,<P) 


T(t-s)(n,<p)  -  |BT(t-O)^((Vs)"1Ws(n>(p))(0)d0. 


(2.4) 


It  can  be  shown  [6]  that  (s(t,s)|0  _<  s,  t  _<  t*}  is  an  evolution  operator  on  Z. 
Taking  a  slightly  different  route  in  the  calculations  after  (2.3)  we  find  that 
u*  is  also  given  by 


u*(t)  -  -C"1^*IKt)S(t)(n/P)  almost  everywhere. 


(2.5) 


* 


6. 


ft* 

with  R(t)(n,<p)  =  T*(t*-t)5*S(t*,t)  (n,<p)  +  T*(o-t)^S(a,t)(n,4>)da,  for 

■'t 

t  €  [0,t*].  Moreover  II(t)  €  Sf(Z,Z)  is  found  to  be  nonnegative,  self-adjoint 
and  is  a  solution  of  the  following  Riccati  integral  equation  in  Z; 


n(t)(n,«p)  =  T*(t*-t)5*T(t*-t)(n,<P)  + 

ft*  1 

+  T*(a-t)[^  -  n(a)  0c  -0*n(a)]T(a-t)(n,<P)do 

Jt 


(2.6) 


for  t  £  [0,t*].  Much  of  the  above  development  is  greatly  facilitated  by 
recent  results  [6]  on  the  regulator  problem  in  a  general  Hilbert  space;  the 
details  for  the  case  of  delay-differential  equations  are  given  in  [9].  Formal 
differentiation  of  (2.6)  leads  to 


jt"W> 

n(t*) 


-A*lKt)z  -  n(t)Az  +  n(t)^C_1-^*Tl(t)z  -  ^z  for 
& 


t  e  [o,t*n 

>(2.7) 


where  z  *  (n,4>)  and  A  and  A*  are  the  infinitesimal  generators  of  the 
semigroups  T(t)  and  T*(t),  respectively.  Of  course,  (2.7)  resembles  the 
well-known  Riccati  differential  equation  arising  in  optimal  control  theory  for 
the  linear  ordinary  differential  equation  x(t)  -  AqX(j;)i  and  many  former 
investigations  have  been  directed  towards  finding  the  "correct"  form  of  the 
analogous  Riccati  equation  for  delay  equations.  This  is  not  an  easy  task, 
since  the  operators  A  and  A*  are  differential  operators  with 

•  m 

Dom(A)  »  (  (<P(0)  ,<P)  |  <P  absolutely  continuous  and  <p  €  L  },  A(q>(0),<p)  ■  (Up,«p); 

for  a  characterization  of  A*  we  refer  to  [ 14] .  For  one  discrete  delay 


■  .  i. 


7. 


(I  =  1,  A_j  =  0)  coupled  systems  of  nonlinear  partial  differential  equations 
of  Riccati-type  in  Euclidean  space  have  been  derived  in  [5,10,12,13].  For 
the  more  general  functional  L  it  is  shown  in  [4]  that  (2.7)  can  be  properly 

studied  in  the  densely  injected  triple  of  spaces  jj^’2  <=.  Z  a  OiT1,2)*, 

12  "2 
where  W  ’  «  {(<P(0),<p)  €  z|  <p  absolutely  continuous  and  cp  €  lr)  and 

I?*  12 

’  )  is  the  dual  of  '  .  The  Riccati  integral  equation  (2.6)  avoids 
many  of  the  technical  difficulties  arising  in  the  development  of  [4]  and  is 
also  more  appropriate  for  our  next  goal,  the  discussion  of  approxima¬ 
tion  schemes  for  (P) .  Finally,  we  mention  that  the  equations  derived  in  this 
section  hold  for  more  general  right-hand  sides  than  those  included  in  L. 

3.  Approximation  Schemes. 

Recent  results  on  the  approximation  of  linear  delay-differential  equations 
allow  discussing  approximation  methods  for  (P)  within  a  general  framework,  including 
a  number  of  examples  of  specific  schemes.  The  objective  is  to  approximate 
(P),  or  equivalently  (^) ,  by  a  sequence  of  optimal  control  problems  for 
ordinary  differential  equations  in  Euclidean  space  (linear  regulator  problems) 
for  which  computer  algorithms  are  readily  available.  We  shall  make  use  of  4 

the  following  definitions  and  hypotheses.  With  {ZN},  N  -  1,2,...  we  denote 
a  sequence  of  linear,  finite  dimensional  subspaces  of  Z,  the  canonical 
projections  from  Z  onto  ZN  are  denoted  by  PN,  and  Q^:  Z  -*•  Z  is  the 
operator  given  by  QQ(n,'0  -  (n,0). 


(HI)  There  exists  a  family  of  semigroups  T1* (t) :  Z  -*■  Z,  for  N  -  1,... 
and  t  €  [0,t*],  such  that 


8. 


(i)  1 1 TN (t )  }|  £  M  exp(u)t)  for  some  M  >  0,  to  £  R  and  t  £  [0,t*] 

(ii)  TN(t)ZN  c  ZN  for  all  t  £  [0,t*] 

(ill)  |T(t)z-T  (t)z|  £  p(N,  z)  for  some  real  valued  mapping  p. 

Of  course,  in  the  examples  that  we  have  in  mind  p  will  tend  to  0 
at  a  certain  rate,  as  N  goes  to 

(H2)  lim  PNz  =  z  for  all  z  £  Z. 

N-x» 

N  n  N 

(H3)  There  exists  a  sequence  of  linear  operators  Q  :  IR  Z  ,  N  =  1,2,...  such 
that 

(i)  ||QN-Qj  „  =  pn(N)  for  some  real  valued  mapping 

0  J^QRn;Z)  Q 
p  with  lim  Pn(N)  =  0, 

^  N-*»  ^ 

HQ1*!!  „  <  q  for  some  q  >  1  independently  of  N. 

^0Rn,Z)  “ 

N 

Conditions  (HI)  and  (H2)  are  satisfied  if  Z  are  chosen  as  the 

subspaces  arising  from  averaging  approximations  [1]  or  spline  functions 

[2];  in  these  cases  lim  p(N,z)  =  0  for  all  z  €  Z.  (H3)  will  be  of 

N 

importance  in  the  approximation  of  (2.2),  which  can  be  written  as 

z(t)  =  T(t)z  +  [  T(t-s)QnBu(s)ds ,  with  z  *  (p,tp).  The  obvious  approach 
Jo  U  ft 

to  approximate  (2.2)  is  to  study  zN(t)  *  TN(t)PNz+  TN(t-s)PNQ  Bu(s)ds 

Jo 

as  N  tends  to  ®.  In  [9]  it  is  shown  that  although  this  leads  to 
satisfactory  convergence  results,  rate  of  convergence  results  for  the 
optimal  control,  payoff  and  trajectory  cannot  be  expected  in  general. 

Thus  we  are  led  to  study  an  approximate  optimal  control  problem  of  a 


more  elaborate  form: 


9. 


(^N,U)<i 


Minimize 
,N+U 


J0(POT^z,u)  =  (PN+yJrPN+VzN,U(t*),z",VJ(t*)  )z  + 


t* 


«PN+V^PN+yzN’M(t),zN*y(t)>_  +  (Cu(t),u(t))  Jdt 


]R 


over  u  €  L^(0,t*;IRra) ,  subject  to 


zN,U(t)  =  TN,y(t)PN+yz  +  [  TN+y(t-cr)QNBu(a)d0,  for  t  6  [0 ,  t*] . 

JO 


It  is  well  known  that  estimates  on  the  approximation  of  a  real  valued  function 

by  piecewise  polynomial  functions  depend  on  powers  of  the  derivative  of  the 

approximated  function  [11];  this  explains,  roughly  speaking,  the  problem  for 

convergence  rates  for  (2.2),  which  contains  the  "jump  operator"  Qq.  In 

(^N,y)  we  try  to  avoid  this  difficulty  by  introducing  the 

N  N 

family  of  operators  Q  .  In  certain  cases  Q  will  act  as  "smoothing 

n  N 

operators"  by  mapping  IR  in  the  subspace  Z  consisting  of  smoother 

N  N 

functions  than  Z.  In  practice  Q  could  be  chosen  to  be  P  Qq,  or  also 

N 

as  some  interpolating  operator  appropriate  for  a  previously  chosen  Z  .  It 

is  quite  clear  that  (j^N»y)  can  be  solved  analytically  in  the  same  fashion 

as  (j^)  and  we  denote  the  optimal  controls  and  optimal  trajectories  of 

(^N,y)  by  GN’y  and  zN,y (t;uN,y) ,  respectively.  In  the  convergence  results 

below  we  shall  use  p(N,z)  ■  p(N,z)  +  M  exp(u>t*)  |  P*z-z  |  for  N  *=  1,2,...  and 

z  €  Z,  and  p(N+h,QN)  =  max  P(N+y,[QN] ,),  where  [QN]  is  the 

j»l,...,N  J 

n  x  n-matrix  representation  of  Q1* :  lRn  ZN,  whose  column  vectors  CQN]  j  are 
elements  of  Z.  Employing  (2.4)-(2.6)  one  can  derive  the  following  convergence 


results: 


10. 


Theorem  3.1.  Assume  that  (H1)-(H3)  hold.  Then  for  each  e  >  0  there  exists 
an  index  and  a  constant  such  that 


JIn,P 
I  u  -u* I 


,  m  <  e  +  K1(p(Nn-HJ,Q  )  +  p(N_+p,z)). 

L2(0,t*jlRra)  1  0  0 


In  the  above  theorem  e  =  K2 (p^(N)+ ||PN+tJQ0-Q0 1|  );  the  constants  and 

K2  ,  and  the  index  Nq  can  be  calculated  explicitly  in  terms  of  the  parameters 
of  (P).  Theorem  3.1  therefore  establishes  that  the  rate  of  convergence  of 
the  approximating  optimal  controls  u^,y  of  (i^^,y)  into  any  specified 
e-neight  orhood  of  u*  is  determined  by  the  rate  of  convergence  of  the  un¬ 
perturbed  semigroups  TN(t)  to  T(t)  when  acting  on  tQN] ^  and  z.  Of 
course,  for  p  =  0  Theorem  3.1  asserts  convergence  of  the  optimal  controls 
UN’°  to  u*  if  only  p(N,z)  and  P(N,[QN]j)  tend  to  0  as  N  goes  to 
infinity. 


Theorem  3.2.  Let  (H1)-(H3)  hold.  Then 


(i)  |z(t;u*)-zN’y(t;GN’y)|  <  K.[p(N+p,z)+p(N+p,QN)+pn(N)+|u*-uN,y|  , 

“  3  Q  L2(0,t*dRm) 


(ii)  |(n(t)z-flN,y(t)PN+yz,y  )z|  <  K4(p(N+P,y) |z|+|y| [p(N+p,^+| z | p(N+p ,QN)  + 


+  |u*-uN,y|  +  |z| 

L2(0,t*-IRra) 


Pq(N)  +  |  z  I  ||Q0-PN+yQ0H  ]), 


for  all  y  and  z  in  Z;  the  constants  and  can  be  calculated  ex¬ 
plicitly  and  does  not  depend  on  y  or  z. 


It  is  simple  to  see  that  Theorems  3.1  and  3.2  imply  an  estimate  on  the 
rate  of  convergence  of  .  JQ(PN+lJz,uN,lJ)  to  JQ(z,u*),  see  [9]. 

Finally,  we  also  have  the  following  result  on  convergence  of  in 

the  supremum  norm: 


Theorem  3.3.  If  (H1)-(H3)  hold  and  if  lim  p(N,z)  =  0  for  all  z  €  Z,  then 

r  N*>c0 

lim  sup  |uN  *  (t)-u*(t)  |  =  0. 

N-x*>  t€  [0,t*] 

To  use  the  above  theorems  for  actual  computations  we  choose  bases 

( )  of  the  H^-dimensional  spaces  Z^.  The  matrix  representations 

of  the  infinitesimal  generators  TN(t)  and  TN(t)*  with  respect  to  the  chosen 

basis  in  Z  are  denoted  by  and  A^,  respectively.  Similarly,  we  use 

[FN]  and  [DN]  for  the  matrix  representation  of  PN^rPN  and  pN.S>pN.  For 

any  pair  of  integers  N  and  y  the  coordinate  vector  of  the  representation 

of  zN’y(t;uN*y)  in  Z^"^1  is  denoted  by  w^’^(t)  and  ir^’^(t)  stands  for 

the  matrix  representation  of  nN»^(t).  Due  to  the  fact  that  QN  maps  ]Rn 
N 

into  Z  we  need  an  additional,  but  not  restrictive,  hypothesis: 


(H4)  For  each  N  ■  1,2,.., 


ZN  c  2** 


there  exists  an  integer  y  >  0  such  that 


For  N  ■  1,2,...  we  assume  that  y  ■  y(N)  is  chosen  according  to  (H4) 
and  we  let  ( QN ] and  [(QN)  ]N+)J  denote  the  matrix  representations  of 
Q*  •  ZN+^  and  (Q^)  :  ZN+^  -*•  ]Rn  with  respect  to  the  bases 

)  of  ZN.  We  are  now  prepared  to  present  the  solution  of 


f 


12. 


n><k^+ 

in  terms  of  an  ordinary  differential  equation  in  IR  y,  coupled  with  a 
matrix-Riccati  differential  equation  and  a  feedback  control  law. 


wN’y(t)  =  AN+ywN’V)  +  [QNlN+yB(t)uN,U(t),  for  0  <  t  <  t*.  ) 

wN’y(0)  =w«>y 

/^(t)  =  -Aj*V.W(t)  -  7TN’y(t)AN+y-[DN+y]  + 

+  TTN'y(t)[QN]N4.vB(t)C(t)B*(t)[((^>*]N+vlTTN’li(t)  \ 


(3.1) 


irN’y(t*)  =  [FN+y] 

uN’y(t)  =  -C_1B*[(QN)*]N+M7rN,y(t)wN’y(t) 


where  wj?,yis  given  by  PW+yz  =  £  g^y(w^,y)..  The  system  (3.1)  is  a  con- 

°  i=l  1 

sequence  of  (2.A)-(2.6),  where  we  replaced  all  operators  by  their  respective 
approximating  operators  and  z  by  P^^z;  since  T^+yZ^,y  c  ZN+y  and  by 
(H3)  and  (HA)  these  equations  are  then  actually  equations  in  the  finite 
dimensional  space  ZN+y. 

If  according  to  Theorem  3.1  we  want  to  determine  Nq  for  a  specified 
e  and  subsequently  take  the  limit  as  y  goes  to  infinity,  then  (HA)  has  to 
be  replaced  by 


,N+y„  = 


S-Hi 


,N+P,N,U, 


(HA*)  for  each  N  =  1,2,...  there  exists  a  nontrivial  sequence  y.  -  y.(N) 
N  N+yk(N) 

such  that  Z  c  Z  ,  for  k  -  1,2,... 

For  subspaces  of  averaging  approximations  and  also  for  subspaces  of 


spline  functions  many  of  the  details  that  are  necessary  to  determine  (3.1)  as  well 
as  estimates  on  the  rate  of  convergence  can  be  found  in  [9]. 
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